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The PHENIX Experiment at the Relativistic Heavy Ion Collider has conducted a survey of fluctu- 
ations in charged hadron multiplicity in Au+Au and Cu+Cu collisions at ^/snn = 22, 62, and 200 
GeV. A universal power law scaling for multiplicity fluctuations expressed as a 2 / \i 2 is observed 
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of the measurement. PHENIX has also measured transverse momentum correlation amplitudes 
in p+p, d+Au, and Au+Au collisions. At low transverse momentum, significant differences in the 
correlations between the baseline p+p and d+Au data and the Au+Au data are presented. 
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1. A Survey of Charged Hadron Multiplicity Fluctuations 

The topic of event-by-event fluctuations of the inclusive charged particle multiplicity in rela- 
tivistic heavy ion collisions has been revived by the observation of non-monotonic behavior in the 
scaled variance as a function of system size at SPS energies [fflj. The scaled variance is defined as 
a 2 /ll, where a 2 represents the variance of the multiplicity distribution in a given centrality bin, and 
jj. is the mean of the distribution. For reference, the scaled variance of a Poisson distribution is 1 .0, 
independent of pL. PHENIX has surveyed the behavior of inclusive charged particle multiplicity 
fluctuations as a function of centrality and transverse momentum in y/sjw = 62 GeV and 200 GeV 
Au+Au and 22, 62, and 200 GeV Cu+Cu collisions in order to investigate if the non-monotonic 
behavior persists at RHIC energies. 

Details about the PHENIX experimental configuration can be found elsewhere [^]. All of the 
measurements described here utilize the PHENIX central arm detectors. The maximum PHENIX 
acceptance of 1 17 1 < 0.35 in pseudorapidity and 180° in azimuthal angle is considered small for 
event-by-event measurements. However, the event-by-event multiplicities are high enough in RHIC 
heavy ion collisions that PHENIX has a competitive sensitivity for the detection of many fluctuation 
signals. For example, a detailed examination of the PHENIX sensitivity to temperature fluctuations 
derived from the measurement of event-by-event mean pj fluctuations is described in [Bt]. 

It has been demonstrated that charged particle multiplicity fluctuation distributions in elemen- 
tary and heavy ion collisions are well described by negative binomial distributions (NBD) [Q]. The 
NBD of an integer m is defined by 

v ; m!(Jfc-l)! (l + il/k) m+k 

where P(m) is normalized for < m < °°, /I =< m >. The NBD contains an additional parameter, 
k, when compared to a Poisson distribution. The NBD becomes a Poisson distribution in the limit 
k^oo. The variance and the mean of the NBD is related to k by l/k = a 2 /jJ. 2 - 1/jU. The PHENIX 
multiplicity distributions are well described by NBD fits for all species, centralities, and transverse 
momentum ranges. The data presented here are results of NBD fits of the multiplicity distributions, 
an example of which is shown in Fig. |]. 

Each 5% wide centrality bin selects a range of impact parameters. This introduces a com- 
ponent to the multiplicity fluctuations that can be attributed to fluctuations in the geometry of the 
collision. It is desireable to estimate and remove this known source of fluctuations so that only fluc- 
tuations due to the dynamics of the collision remain. The contribution of geometrical fluctuations 
is estimated using the HIJING event generator [^j], which well reproduces the mean multiplicity of 
RHIC collisions [g]. The estimate is performed by comparing fluctuations from events with a fixed 
impact parameter to fluctuations from events with a range of impact parameters covering the width 
of each centrality bin, as determined from Glauber model simulations. The HIJING estimates are 
confirmed by comparing the HIJING fixed/ranged fluctuation ratios to measured l%/5% bin width 
fluctuation ratios. The magnitude of the scaled variance is reduced by this correction. A 15% 
systematic error for this estimate is included in the errors shown. 

By measuring the scaled variance in successively wider azimuthal ranges, a linear dependence 
on azimuthal acceptance is observed, with the scaled variance increasing with azimuthal accep- 
tance. In order to facilitate direct comparisons with other experiments, the multiplicity fluctuations 
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Figure 1: PHENIX Preliminary. The event-by-event inclusive charged hadron multiplicity distribution for 
0-5% central 200 GeV Au+Au collisions with transverse momentum in the range 0.2 < pj < 2.0 GeV/c. 
The dashed line is a Negative Binomial Distribution fit to the data. 



quoted here have been linearly extrapolated to 2n acceptance by fitting the azimuthal dependance 
within the detector acceptance. Systematic errors due to the extrapolation have been included in 
the total errors shown. The results of the survey are shown in Fig. || through ||, where the corrected 
scaled variance is plotted as a function of centrality and transverse momentum range. 

In the Grand Canonical Ensemble, the variance and the mean of the particle number can be 
directly related to the compressibility: d 2 //i 2 = ks(T /V)kj , where ks is Boltzmann's constant, 
T is the system temperature, and V is its volume [j7j]. The multiplicity fluctuations in terms of 
a 2 / \l 2 = 1 /k^BD + 1 /jU are shown in Fig. ^ as a function of centrality for Au+Au collisions at 200 
and 62 GeV and Cu+Cu collisions at 200, 62, and 22 GeV. In order to demonstrate their scaling 
properties as a function of centrality, each curve has been scaled to the 200 GeV Au+Au curve. All 
four datasets can be described by a power-law function as follows: a 2 / ' jl 2 « AC,^ . 

It is expected that the compressibility diverges as one approaches the critical point, and the rate 
of divergence is described by a power law, k T = A((T — T c ) / T c )~ r , where Tq is the value of the 
temperature at the critical point, and y is the critical exponent for isothermal compressibility [Q]. 
However, it is not prudent to interpret the observed scaling as critical behavior since it has not been 
established that the temperature varies significantly as a function of N part . A more robust critical 
exponent study could be performed with a RHIC scan at low energy where the freeze-out tempera- 
ture could be estimated from identified particle spectra and Hanbury-Brown Twiss correlations, as 
outlined in [0]. 

Figures ^ and [7] show the multiplicity fluctuations as a function of centrality for different 
Pt ranges. In each case, the exponent of the power law function describing the data remains 
unchanged. Hence, the influence of pj -dependent processes at higher pj, such as hard scattering, 
have little effect on the scaling properties. The scaling also applies to the low energy 22 GeV 
Cu+Cu data. 
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2. Low Transverse Momentum Correlations 

The technique of two-particle correlations has demonstrated its effectiveness in extracting in- 
formation about jet production in reference (p+p and d+Au) and heavy-ion collisions at RHIC. 
Some features of the correlations are dependent on the pj of the trigger and associated particles. 
It is important to examine and understand the evolution of correlation functions from high to low 
Pt in heavy ion collisions with comparisons to reference collisions. Examining the pj evolution 
of the correlations may provide information that can aid in the interpretation of such features such 
as excess jet production at low pj or eventually the observed displaced away side peaks. 

Interpreting correlation functions at low pj (in this context, with both the trigger and associ- 
ated particle below 1 GeV/c) is complicated by several factors: a) the hard scattering correlation 
becomes reduced in amplitude and diluted in azimuthal separation, b) correlations due to collective 
flow are present, although the magnitude of elliptic flow is much reduced at low pj, c) Hanbury- 
Brown Twiss correlations are also present, and d) correlations due to resonance production are 
present. In order to help disentangle the latter two contributions, all correlations are separated into 
only like-sign or only unlike-sign pairs. Here, only the like-sign results, which should not have a 
significant resonance contribution, are presented. 

The minimum bias datasets analyzed here contain 10 million events from the RHIC Run-3 p+p 
run, 10 million events from the Run-3 d+Au run, 2 million events from the Run-4 200 GeV Au+Au 
run, and 750,000 events from the Run-4 62 GeV Au+Au run. All correlation functions are con- 
structed using a strict mixed event procedure with all cuts, including track proximity cuts, applied 
identically to the data and mixed event samples. Mixed events are constructed from identical data 
classes (distinguished by centralities within 5% and event vertex z coordinates within 10 cm) by 
directly sampling the data charged particle multiplicity distribution after all cuts have been applied. 

Since every event contains multiple particles at low pj, the correlations are constructed by 
considering every particle in the event as a trigger particle in succession and associating it with 
every other particle in the event. Correlations constructed in this manner are auto-correlations [|10|]. 
The correlation amplitudes reported here are normalized bin-by-bin as follows: C = ® iau <l^" ems - dau < t 

™ mixed I ™ events, mixed 

where Ndata and N m i xe d are the number of pair counts in a given pj and/or A(j) bin. Correlations 
shown for the near-side are defined as the integrated correlation amplitude over the range |A0| < 
60°. 

Correlation amplitudes plotted as a function of the transverse rapidity, yj = ln{{mj +pr)/ino) 
with the pion mass assumed for rriQ, for like-sign pairs on the near-side, are shown in Fig. [8] for 
minimum bias 200 GeV d+Au collisions and in Fig. || for 0-5% central 200 GeV Au+Au collisions. 
The yj variable is chosen in order to emphasize the distribution at low pj with respect to that at high 



pj and to facilitate comparisons to measurements by the STAR Collaboration (jlOJ. For reference, 
yr = 1.5 corresponds to pj = 300 MeV/c and yj = 2.7 corresponds to pj = 1.0 GeV/c for pions. 
There are two primary features observed in these correlations. First, there is the expected large 
correlation amplitude primarily due to hard scattering processes when either particle in the pair 
has a large pj along with a small fractional contribution due to elliptic flow. Second, there is a 
peak when both pairs have low pj. This peak is dominated by HBT correlations, which has been 
confirmed by observing a sharp reduction in the amplitude of the peak for unlike-sign pairs and 
by simultaneously observing a peak in the Qi m ariam distribution in this region for like-sign pairs. 
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Closer xamination of the low pj peak reveals that there is a difference in the location of the peak 
in yj between d+Au and Au+Au collisions. In d+Au collisions, the peak amplitude rises up to the 
low pt edge of the PHENIX acceptance. This holds for p+p collisions, also. However, there is a 
maximum in the peak at yj 1 .4 (pj ~ 250 MeV/c) in Au+Au collisions at both 200 and 62 GeV. 
Fig. 10 shows the peak taken along the pj \ = pj2 diagonal for near-side like-sign pairs. Possible 
explanations for these differences require further study. 

To facilitate a more detailed direct comparison of the d+Au and Au+Au correlations over the 
entire pj range, it is necessary to subtract the flow component from the Au+Au azimuthal corre- 
lation distributions. This is done by fitting the near-side region (A0 < 90°) of these distributions 
with a function containing a harmonic cos(2A(f)) term and a Gaussian distribution with a mean of 
Aip = 0. For p+p and d+Au collisions, the Gaussian fit is performed, but the harmonic term is not 
included. 

Fig. [n] shows the correlation amplitude from the near-side Gaussian fit for like-sign pairs as a 
function of centrality. For this plot, all pairs have a pj between 200 and 500 MeV/c. In order to best 
isolate the HBT contribution to the peak, all pairs are required to be within |Atj| < 0.1. For both 
200 and 62 GeV Au+Au collisions, the amplitude drops exponentially as a function N part i C ip ants . 
Only a small fraction of the decrease in the amplitude is due to random dilution of the correlation 
amplitude [JTXjj. 

Fig. |l2| shows the standard deviation of the near-side Gaussian fit for like-sign pairs as a 
function of centrality. Again, both particles in all pairs have a pj between 200 and 500 MeV/c and 
|Atj| < 0.1. The error bars are dominated by systematic errors due to the fitting procedure. There 
is a large difference observed in the width between p+p and d+Au collisions. In Au+Au collisions, 
particularly for the higher statistics 200 GeV Au+Au dataset, a significant increase in the width is 
seen in the most central collisions. 



Fig. 13 shows the correlation amplitude for the near-side Gaussian fit for like-sign pairs as 
a function of transverse momentum. Here, pT,min < Pt,i < Pr,max and pr,min < Pt,2 < Pr,max (Pt 
within the limits of the bin, pT,mtn and pr,max) and |Atj | < 0.1. The points are plotted in the geomet- 
ric center of the bin. The bins are exponentially distributed in order to nullify any variations due 
to random dilution of the correlation amplitudes within each collision species. In both the p+p and 
d+Au datasets, the amplitudes are relatively flat as a function of pj up to about 1-1.5 GeV/c, where 
the amplitudes begin to rise significantly. This rise is likely due to the influence of correlations due 
to jets, which has been confirmed by the appearance of an away-side component in the azimuthal 
correlations. A rise in the correlation amplitudes for Au+Au collisions is also seen starting at the 
same pj. However, unlike the p+p and d+Au datasets, the Au+Au data exhibit an exponential 
decrease in the amplitude as a function of pj below 1 GeV/c. This decrease is contrary to the trend 
in the reference datasets, and contrary to the ansatz of significant contributions from jet production 
at low pj [10]. More detailed studies of the cause of this effect are underway. 



3. Summary 

PHENIX has completed a comprehensive survey of charged hadron multiplicity fluctuations 
for several different collision species and energies, all of which demonstrate a universal power-law 
scaling as a function of N part . Also presented is a survey of two-particle correlations for several 
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collision species at 62 and 200 GeV focusing on the behavior of the near-side peak for like-sign 
pairs, which is dominated by HBT correlations. Examination of spatial HBT correlations shows 
several statistically significant differences between p+p and d+Au collisions and Au+Au collisions 
that will warrant further study. 
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Figure 2: PHENIX Preliminary multiplicity fluctua- 
tions for inclusive charged hadrons in terms of <y 2 / jX 
as a function of N 'participants and the transverse mo- 
mentum range 0.2 < pj < pr.max GeV/c for 200 GeV 
Au+Au collisions. The data have been corrected to 
remove impact parameter fluctuations and have been 
extrapolated to 2n azimuthal acceptance. 




Figure 4: PHENIX preliminary multiplicity fluctua- 
tions for inclusive charged hadrons in terms of a 2 j[l 
as a function of N participants and the transverse mo- 
mentum range 0.2 < pj < pr,max 

GeV/c for 200 GeV 
Cu+Cu collisions. The data have been corrected to 
remove impact parameter fluctuations and have been 
extrapolated to 2k azimuthal acceptance. 




Figure 3: PHENIX preliminary multiplicity fluctua- 
tions for inclusive charged hadrons in terms of <J 2 /n 
as a function of N 'participants and the transverse mo- 
mentum range 0.2 < pj < pr.max GeV/c for 62 GeV 
Au+Au collisions. The data have been corrected to 
remove impact parameter fluctuations and have been 
extrapolated to 2% azimuthal acceptance. 




Figure 5: PHENIX preliminary multiplicity fluctua- 
tions for inclusive charged hadrons in terms of a 2 / ji 
as a function of N 'participants and the transverse mo- 
mentum range 0.2 < pj < pr.max GeV/c for 62 GeV 
Cu+Cu collisions. The data have been corrected to 
remove impact parameter fluctuations and have been 
extrapolated to 2n azimuthal acceptance. 
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Figure 6: Multiplicity fluctuations for inclusive 
charged hadrons in the transverse momentum range 
0.2 < pt < 2.0 GeV/c in terms of a 2 / [l 2 as a function 
of Nparticipants f° r 200 GeV Au+Au, 62 GeV Au+Au, 
200 GeV Cu+Cu, 62 GeV, and 22 GeV Cu+Cu col- 
lisions. The data have been scaled by factors of 1.0, 
0.73, 1.43, 0.73, and 1.0, respectively, in order to 
emphasize the universal scaling of all species. The 
dashed curve is a power law fit as described in the text. 



Figure 7: Multiplicity fluctuations for inclusive 
charged hadrons in the transverse momentum range 
0.2 < p T < 0.75 GeV/c in terms of G 2 //i 2 as a func- 
tion of Np ar ticipunts for 200 GeV Au+Au, 62 GeV 
Au+Au, 200 GeV Cu+Cu, 62 GeV, and 22 GeV 
Cu+Cu collisions. The data have been scaled by fac- 
tors of 1.0, 0.73, 1.43, 0.73, and 1.0, respectively, in 
order to emphasize the universal scaling of all species. 
The dashed curve is a power law fit as described in the 
text. 



PHENIX Preliminary: 200 GeV d+Au, Min. Bias, Like-Sign, Near-Side 




Figure 8: Near-side correlation amplitudes for like- 
sign pairs as a function of the transverse rapidity of 
each particle in the pair for minimum bias 200 GeV 
d+Au collisions. 



PHENIX Preliminary: 200 GeV Au4.Au, 0-5% Gentrality, Like-Sign, Near-Side 




Figure 9: Near-side correlation amplitudes for like- 
sign pairs as a function of the transverse rapidity of 
each particle in the pair for 0-5% central 200 GeV 
Au+Au collisions. 
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Figure 10: Near-side correlation amplitudes for like- 
sign pairs with pj bins chosen such that pj\ — Pt,2 

for 

Au+Au and d+Au collisions. These correlation ampli- 
tudes include contributions from flow in the Au+Au 
collisions. The error bars include statistical and sys- 
tematic errors. 
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Figure 11: Near-side correlation amplitudes for like- 
sign pairs with pj between 200 and 500 MeV/c and 
\Ari\ < 0.1 for p+p, d+Au, and Au+Au collisions. 
Flow has been subtracted from the Au+Au data. The 
errors are dominated by systematic errors due to the fit 
procedure. 
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PHENIX Preliminary, Like-sign pairs 
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PHENIX Preliminary. Like-sign pairs. 
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Figure 12: Near-side correlation widths for like-sign 
pairs with pj between 200 and 500 MeV/c and \Arj | < 
0.1 for p+p, d+Au, and Au+Au collisions. Flow has 
been subtracted from the Au+Au data. The errors are 
dominated by systematic errors due to the fit proce- 
dure. 



Figure 13: Near-side correlation amplitudes for like- 
sign pairs both with pr within the given pr bin and 
with |At}| < 0.1 for p+p, d+Au, and 0-5% central 
Au+Au collisions. Flow has been subtracted from the 
Au+Au data. The errors are dominated by systematic 
errors due to the fit procedure. 
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